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We study the fractional quantum Hall effect at i^ = 7/3 and 5/2 in the presence of the spin-orbit 
interaction, using the exact diagonalization method and the density matrix renormalization group 
(DMRG) method in a spherical geometry. Trial wave functions at these fillings are the Laughlin state 
and the Moore-Reed-PfafBan state. The ground state excitation energy gaps and pair-correlation 
functions at fractional filling factor v — 7/3 and 5/2 in the second Landau level are calculated. We 
find that the spin-orbit interaction stabilizes the fractional quantum Hall states. 
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The fractional quantum Hall (FQH) effect is a re- 
markable many-body phenomenon observed in two di- 
mensional electrorL-gases (2DEG) in a strong perpendic- 
ular magnetic fielcLl. The Laughlin state and it's hierar- 
chy theory describe the fractional quantum Hall (FQ|HQ. 
effect at odd-denominator filling factors successfullyBtl 
The feature of Laughlin state is the appearance of quasi- 
particle excitations with fractional charge and fractional 
statistics called anyon.cl On the other hand, the Laughlin 
state does not include an even-denominator state. There- 
fore the quantum Hall plateau observed at filling factor 
1/ — 5/2 (second Landau level even dominator) poses a 
special challenge. 

There are some candidates for the ground state (GS) 
aX V = 5/2. Particularly, the Moore-Read Pfaffian (Pf) 
state is believed to be a strong caxididate for fractional 
quantum Hall state at this fillingJj'tl The Pfaffian wave 
function is written by 
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iyi, the Pfaffian factor is defined by 

^ EaGS« Sg^^ nfi? A//<T(2fc-l)a(2fc) , 



^ij ~ 2«/2(7V/2)! 

for an N X N antisymmetric matrix whose elements are 
Mij-, Sn is the group of permutations of A^ objects. The 
Pfaffian state is totally antisymmetric for even denomi- 
nator, so could describe electrons.Q 

The occurence of the Pf state is very exciting since 
the quasiparticle excitations in this state has non-abelian 
statistics (Laughlin state excitations are abelian anyons) . 
The non-abelian topological phase suggests potential ap- 
plications toward quantum computing. The quantum 
gates constructed from abehan anyons are very limitedJj 
Non- Abelian anyons, on the contrary, are much more 
useful to the topological quantum computer,Ll since the 
braiding of non-Abelian anyons induces noncommuting 
processes, one can construct various quantum gates from 
non-Abelian representations of braid group. Therefore, 
determining the nature of i^ = 5/2 state is of greater 
urgency, beyond the conventional FQH physics. 



Unfortunately, experiments showed that the energy 
gap at the v = 5/2 is vpj^mall and then the FQH state 
at this filling is fragileEHlEJ. Numerical studies indicate 
the fact that the Pf state exist between the compress- 
ible stripe state and the corrmnsit fermion liquid (CFL) 
state, in the phase diagram.E£rEj These studies showed 
the v = 5/2 ground state is near of the stripe phase and, 
extra short range interaction changes Pf state into stripe 
state. 

In this paper we show that the Pf state could be more 
stable by moving the point away from the transition 
point, and the spin-orbit interaction has good influence 
on stability of the pfaffian state. A large number of theo- 
retical studies have been devoted to investigate the FQH 
states in 2DEG, however, the effect of the spin-orbit in- 
teraction on the FQH state has not been studied enough 
yet. In this paper we investigate the FQH state in the 
presence of the spin-orbit interaction by using numerical 
method. In 2DEGs, the spin-orbit interaction has two 
relevant coplpibutions, the Rashba term and the Dressel- 
haus term.By 

The Rashba coupling Hamiltonian is described by 

HB.^x{n,ay-nya^). (2) 

The Dresselhaus Hamiltonians is described by 



HD = X'{Il.,a'=^Uyay). 
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Here, the spin-orbit coupling parameter A is linearly de- 
pendent on the expectation value of the interfacet^ectric 
field (Ez) at 2DEG, this relation is described byO 



A^(i?.), 



(4) 



here (Ez) is controlled by the gate voltage. Rashbha 
spin-orbit coupling arises from the structure inversion 
asymmetry of the quantum well, and the Dresselhause 
term stems from the bulk-inversion asymmetry of semi- 
conductor material. The coefficient A' of the Dressel- 
hause term is fully determined by the geometry of the 
heterostructure, while the Rashba coefficient A can be 
varied by an electric field across the well. Lommer ct al. 



have theoretically pointed out that the Rashba mecha- 
nism becomes dominant in a narrow gap semiconductor 
system.ta Later, Luo et al. have experimentally shown 
that the Rhashba term is donuaant for spin splitting in 
an InAs based heterostructureJlj We take account of only 
the Rashba term in this paper. 

Eigenstates of the single-particle Hamiltonian including 
the Rashba term can be written in the form 

1-0) — cos (—a) \n = 1) + sin ( —aj \n = 0), (5) 

where \n = 1) is a second Landau level (SLL) eigenfunc- 
tion, \n = 0) is a lowest Landau Level (LLL) eigenfunc- 
tion, cos(a) is written byj^ 



"K \ /I hWc/4: 
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From eqs. (4) and (6), we find that the mixing parameter 
a increases from a = with the increase of the perpen- 
dicular electron field. In the limit of strong spin-orbit 
coupling or weak magnetic field, the value of a is 0.5 



(cos(|a) 



This limit corresponds to the n = 1 



LL of graphcncEil We note that since the GS at v ^\/2 
in the LLL(a = 1.0) is shown to be the CFL statc,El the 
Pf state will be stabilized by the weak spin-orbit interac- 
tion. 

To study many-body problems, we study the projected 
many body interaction Hamiltonian onto a certain Lan- 
dau Level which is written as,u 
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where Ri is the guiding center coordinate of the 
ith electron, which satisfies the commutation relation, 
[i?J,-R^] — ii%Sjk, and V{q) = '^ttj- is the Fourier trans- 
form of the Coulomb interaction. In this work, the mag- 
netic length Ib is set to be 1, and we take e^/eis as units 
of energy scale. We omit the component at g = 0, which 
is canceled by uniform positive background charge. We 
assume that the width of the wave function perpendicu- 
lar to the two dimensional plane is-fiufficiently small com- 
pared with the magnetic length.E3 The relativistic form 
factor in the nth LL is written as 



Fo{q) = Lo(y), 



(8) 



and 



F„>i(q) =cos2(-a)L„(|-) + sin2(-a)i„_i(:|-), (9) 

here, L„(x) are the Laguerre polynomials. 

In this paper, we investigate the effects of the spin- 
orbit interaction on the Laughlin state and the Pfaf- 
fian state by applying this projected Hamiltonian. We 
calculate the exact wave functions of the ground state 



for several jiialujes of a, using the exact diagonalization 
for N < 14,00 and the density matpix renormalization 
group (DMRG) method for N <18E3E3, in the spherical 
geometry. The DMRG method is a real space renormal- 
ization group method combined with a exact diagonal- 
ization method. The DMRG method provides low-energy 
eigenvalues and corresponding eigenvectors of Hamiltoni- 
ans within a restricted number basis states. The accuracy 
of the results is systematically controlled by the trunca- 
tion error, which is smaller than 0.0003 in the present 
calculation. 

In the spherical geometry, it is convenient to write the 
Hamiltonian aaJ 
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(10) 



where P^ [771] projects onto states in which particles i and 
J have the relative angular momentum hm, and Vm{n) is 
their interaction energy in the nth LL. Using the above 
form factors, the pseudopotentials arc given by 
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To extrapolate the energy gaps, we take the finite size 
correction to V^ for each system size.Es 

We first report the results obtained by the exact 
diagonalization method at i' = 7/3.CJ Although the 
maximum value of a given by the spin orbit interaction is 
0.5, we increase a for a > 0.5 to see the relation between 
the Laughlin state and the mixed state written by eq. 
(5). Figure 1 shows the energy gaps for 0.2 ^ a ^ 1.0.^^^ 
At a = 1.0, the value of the extrapolated gap obtained 
by the best linear fit is 0.101. rThis value is equal to 
the result of the previous work.t3 The energy gap first 
increases and then decreases with the increase of a. The 
a which gives the maximum energy gap is around 0.55. 
We discuss this behavior by considering the relation 
between the pseudopotential and the energy gap in the 
Laughlin state. The energy gap of the Laughlin state 
is roughly proportional to the difference of V^ and V^ 
obtained by eq. (11). V^ — V^ is written as a function 
of a by 



Vi-y^: = -0.25963 cos''(^a)-|-0. 19386 cos2(^a)+0. 16616. 

(12) 
We compare V^ — V^ with the value of the extrapolated 
energy gaps in Fig. 2. This figure shows that the a de- 
pendence of the energy gap and the a dependence of 
Vi — V^ are similar. This result also means that the GS 
for 0.2 ^ a ^ 1.0 is expected to be the Laughlin state. 
Therefor the effect of the spin-orbit interaction on the 
Laughlin state is featured by the change in the pseudo 
potential. We therefor conclude that the Laughlin state 
is stabilized by the spin-orbit interaction. 

Next we investigate the GS ai v = 5/2 to see the effect 
of the spin-orbit interaction on the Pf state. In Fig. 3, 
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FIG. 1: Energy gaps a,t u — 7/3. The straight lines denote the 
best linear (in 1/A'^) fit to the data points. N is the number 
of electrons. 
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FIG. 2: The cross marks represent extrapolated energygaps. 
The broken curve denotes Vi — V3 . The horizontal axis is 
cos(-|a). 



wc show our results of the energy gaps for ^ a ^ 0.3. 
Half the sum of quasiparticle and quasihole excitation 
energies are plotted as a function of 1/N, where N is 
number of electrons in the system. At a = 0, the plotted 
data are well fitted by a straight line in 1/A^. Our result 
at a = calculated by the DMRG method is 0.025. This 
value is equal to the prewous result calculated by exact 
diagonalization methodJ23 The results in Fig. 3 show that 
the energy gap first increases and then decreases with the 



increase of a from a ~ 0. 

We discuss the a dependence of the energy gap by 
considering the phase diagram of the GS. In general, 
the value of the energy gap depends on the distanpa 
from the phase boundary. In the previous studies, E£ro 
an extra increase of the pseudopotential induce the 
transition from the PfafSan state to the CFL phase. 
The increase of the pseudopotential corresponds to the 
increase of a for ^ a ^ 0.3. The GS at a = 0.1 is away 
from the stripe phase and the CFL phase, consequently 
the energy gap at a = 0.1 is larger than the gap at a = 0. 
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FIG. 3: Energy gaps at p — 5/2. The straight lines denote the 
best linear (in 1/A^) fit to the data points. N is the number 
of electrons. 



In order to further study of the GS around a - 
we next investigate the pair-correlation functions, 
pair-correlation function is defined by 
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Figure 4 shows the results of the pair-correlation func- 
tions for N =18. The curves for a = , 0.1 and 0.2 are 
almost indistinguishable. However, the correlation func- 
tion at a = 0.3 looks that of the CFL state. Precisely 
in a short distance (r < 2), the correlation functions for 
a < 0.2 have shoulder structure which is a signature of 
the Pf state but the correlation function at a = 0.3 does 
not have such a Pf-like structure. In a long distance 
(r > 4), the amplitude of oscillation of the correlation 
function at a = 0.3 is smaller than other correlation func- 
tions. From these differences, we expect that the increase 
in a would induce the phase transition. This result agree 
with our expectation. 

To confirm the Pf-like state around a = 0.1, we finally 
calculate the overlap between the exact ground state and 



the Pf wave function for TV = 6. The PfafBan wave func- 
tion on the spherical geometry is written by 
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here, u = cos(|)e*'^, and v = sin(|)e^*'^, and we choose 
a gauge such that A = (■^) cot(0)</> where 6*, are the 
polar coordinates on the unit sphere with A^$ units 
of flux passing through the sphere. The value of the 
overlap is defined by |(*Gxact|^Pf)|- The overlap at 
a = is 0.866 and it increases with the increase in a. 
The maximum value of the overlap is 0.901 at a = 0.15. 
This behavior is consistent with that of the energy 
gap. We therefor confirmed that the GS at a = 0.1 is 
characterized by the Pf state which is stabilized by the 
spin-orbit interaction. 

In summary, we have investigated the effect of the 
spin-orbit interaction on the Laughlin state and the Pf 
state. We have introduced the parameter a which is in 
the absence of the spin-orbit interaction and it increases 
with the increase of the spin-orbit interaction. For the 
Laughlin state at :^ = 7/3, we find that the energy gap 
is increased by the spin-orbit interaction. The enhance- 
ment of the energy gap is reasonably understood by the 
Haldene's pseudopotential written in eq. (11). For the 
Pf state ai V ~ 5/2, we have calculated the energy gap, 
the correlation functions, and the overlap as a function 
of a. Both the extrapolated energy gap and the overlap 
first increase and then decrease with the increase of a 
from a = 0. The energy gap has the maximum value 
at around a = 0.1 and the maximum of the overlap 
is 0.901 at a = 0.15 for iV = 6. The pair-correlation 
functions for ^ a ^ 0.2 are almost the same. These 
correlation functions are characterized by the Pf-like 




FIG. 4: The pair correlation functions at a = 0, 0.1, 0.2 and 
0.3 for iV = 18. 

state, and we can expect that the GS is the Pf state. 
The pair-correlation functions at a = 0.3 is CFL-like 
state. This means that the increase in a induce the 
phase transition. From above results, we conclude that 
the Pf state is stabilized by the increase of the spin-orbit 
interaction. Spin-orbit coupling linery depends on the 
interface-jelectron field which is controlled by the gate 
voltage.llZl This means that the FQH state is stabilized 
by the gate voltage which is applied perpendicular to 
the quantum well. 
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